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ASSESSMENT OF SITE RESPONSE 


ANALYSIS PROCEDURES 
by 
I. M. Idriss 
1.0 INTRODUCTION 


The Loma Prieta earthquake occurred on October 17, 1989 at 5:04 pm Pacific daylight 
time along a 45-km long segment of the San Andreas fault in the Santa Cruz Mountains. 
The earthquake was assigned a surface wave magnitude, M, = 7.1, and a moment 
magnitude, M,, = 7 making it the largest earthquake to occur in Northern California since 
the 1906 earthquake. 


The earthquake triggered by far the largest number of instruments ever triggered by an 
earthquake and recordings were obtained at well over 200 locations, including free-field 
stations, small buildings, high rise structures and dams. Of these, 88 stations were at 
locations that can be considered to represent free-field conditions, i.e., the instrument had 
been placed in a small instrument shack or is at the ground floor of a small structure 
(three stories or less) with no basement. 


The records obtained from the free-field stations are summarized in Appendix A of this 
report. These include 31 stations at rock sites, nine stations at soft soil sites and 48 
stations at other soils sites. The recordings at the soft soil sites are the largest set of 
recordings ever obtained at such site conditions at these levels of shaking. 


The records obtained at the soft soil sites offered an excellent opportunity to assess the 
procedures used for calculating the response of soil sites during earthquakes. This was 
particularly the case because rock outcrop motions were also available from nearby 
locations and subsurface information was gathered at most of these soft soil sites. In 
particular, six soft soil sites and the rock site at Yerba Buena Island were investigated by 
the Electric Power Research Institute (EPRI) with technical support from USGS and 
partial financial support from the University of California at Davis (through a grant from 
the four companies listed in Section 6.0) and from the San Francisco International Airport 
Authority. The investigation at each soft soil site consisted of drilling through the soil 
profile and into bedrock, some sampling and measurement of shear wave and 
compression wave velocities 


The results of these subsurface investigations at Treasure Island site and at the San 
Francisco Airport site were used in this study to calculate the horizontal components of 
site response at these two sites and to assess the procedures used for conducting such 
response calculations. 


The locations of these two sites and other recording stations pertinent to this study are 
shown in Fig. 1-1. The three-digit numbers shown in the figure identify the recording 
stations operated by the California Strong Motion Instrumentation Program (CSMIP); the 
stations pertinent to this study are the following: 


Station Number Station Name Site Condition 
G7 Treasure Island Soft Soil Site 
223 San Francisco Airport Soft Soil Site 
539 Sierra Point Rock Site 
130 Diamond Heights Rock Site 
131 Pacific Heights Rock Site 
162 Cliff House Rock Site 
133 Telegraph Hill Rock Site 
Tel. Rincon Hill Rock Site 
163 Yerba Buena Island Rock Site 
222 Presidio Rock Site 
not shown in Fig. 1-1 Golden Gate Rock Site 
not shown in Fig. 1-1 Piedmont Rock Site 
not shown in Fig. 1-1 Lawrence Berkeley Lab Rock Site 


The horizontal components of motion recorded at Yerba Buena Island (Station No. 163 in 
Fig. 1-1) were used as rock outcrop motion in the analysis of the horizontal response of 
the Treasure Island site (Station No. 117 in Fig. 1-1) and those recorded at Sierra Point 
(Station No. 539 in Fig. 1-1) were used as rock outcrop motion in the analysis of the 
horizontal response of the San Francisco Airport site (Station No. 223 in Fig. 1-1). The 
horizontal components recorded at the other rock sites listed above were also used as rock 
outcrop in the response calculations for both sites. The locations of the three stations not 
shown in Fig. 1-1 are as follows: the Golden Gate station is just northwest of the Presidio 
station; the Piedmont station is about 7 km northeast and the Lawrence Berkeley 
Laboratory is about 14 km north of downtown Oakland, respectively. 


The characteristics of the horizontal components of the recordings obtained at these rock 
sites are summarized in Appendix B and those obtained at soft soil sites, which include 
the sites at Treasure Island and at the San Francisco Airport and seven other stations, are 
summarized in Appendix C. 


The results of the analyses completed for the Treasure Island and the San Francisco 
Airport are presented in the following sections of this report. 
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Fig. 1-1 Locations of CSMIP Recording Stations in the San Francisco- 
Bay Area used in Present Studies 
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2.0 SITE RESPONSE ANALYSES 


Procedures for calculating the response of soil profile to earthquake shaking have been 
available for a long time. The currently available procedures use either an equivalent 
linear or a nonlinear representation for the stress-strain behavior of the various soil layers 
comprising the soil profile. While other formulations are available, many situations are 
reasonably approximated by a horizontal multi-layered system of semi-infinite extent. 
The latter is most conveniently represented by a one dimensional system consisting of a 
series of soil sublayers over a half-space. 


The equivalent linear procedure was introduced by Idriss and Seed (1968) and has been 
implemented in several computer programs. Currently, the most widely used computer 
program for a multi-layered system over a half-space, which also incorporates the 
equivalent linear formulation, is program SHAKE (Schnabel et al, 1972). Another 
computer program that is also used in many applications is program RASCALS; the 
original program (which was designated RASCAL) was developed and coded by Silva 
and Lee (1987) has been considerably updated by Silva (personal communication, 1992). 


The program SHAKE91 (Idriss and Sun, 1992) was used in the present study. The 
program SHAKE91 is a modification and an update of the original SHAKE program 
(Schnabel et al, 1972). 


There are many nonlinear procedures available for calculating the response of a multi- 
layered system overlying a half-space. For the present study, only two of these 
procedures were examined. One procedure is that incorporated in the program DESRA-2 
(Lee and Finn, 1978) which was developed by Finn and his colleagues at the University 
of British Columbia. The other procedure used in this study is that incorporated in the 
program SUMDES (Li et al, 1992) which was developed by Li and his colleagues at the 
University of California at Davis. 


As part of this study, detailed response analyses were conducted for the Treasure Island 
Site and for the San Francisco Airport Site. Equivalent linear as well as nonlinear 
analyses were completed for the Treasure Island Site and only equivalent linear analyses 
were completed for the San Francisco Airport Site. The results are presented in Sections 
3.0 and 4.0 of this report. 
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3.0 TREASURE ISLAND SITE 


The location of the recording station at Treasure Island is shown in Fig. 1-1 and in more 
detail in Fig. 3-1. Also shown in Fig. 3-1 is the location of the recording station at Yerba 
Buena Island which is the closest rock site to the Treasure Island recording station. The 
spectra of the horizontal motions recorded at Treasure Island and at Yerba Buena Island 
are presented in Fig. 3-2. The upper part of Fig. 3-2 shows the spectra for the stronger 
components (EW components) of the recorded motions and the lower part of the figure 
shows the spectra for the weaker components (NS components). More details regarding 
the characteristics of the recordings at Yerba Buena Island are given in Appendix B and 
those for Treasure Island are summarized in Appendix C. 


3.1 Subsurface Conditions at The Treasure Island Site 


The subsurface conditions at the Treasure Island site were obtained by drilling a boring 
nearby the station to a depth of 338 feet (103 meters). The shear wave velocities were 
measured in the same boring by USGS (Gibbs et al, 1992) and by Redpath (1991) using 
down-hole measuring techniques and by Rollins et al (1992) at other, but nearby, 
locations using a seismic cone (SCPT). 


The log of boring and the shear wave velocities measured by USGS are shown in Fig. 3-3 
and the shear wave velocities measured by USGS, by Redpath (both the originally 
reported as well the revised velocities) and by Rollins et al are shown in Fig. 3-4. Also 
shown in Fig. 3-4 are the values of the “best estimate” shear wave velocities used in the 
analyses; these values are listed below in Section 3.2. 


The log of boring indicates that the subsurface conditions at the Treasure Island site 
consist of sandy fill overlying a natural deposit of loose sand to a depth of about 40 ft 
(12.2 m), which in turn is underlain by a 55-ft (16.8 m) layer of Young Bay Mud. 
Alternating layers of dense sand and Old Bay Mud are encountered below a depth of 95 ft 
(29 m) to a depth of about 290 ft (88.4 m). Weathered shale extends from this depth to a 
depth of 320 ft (97.5 m), where a more competent sandstone is encountered. 


3.2 Dynamic Soil Properties 


The dynamic soil properties needed for an equivalent linear analysis include the 
maximum shear wave velocity (or maximum shear modulus), modulus-reduction as a 
function of strain and damping ratio as a function of strain. The maximum shear wave 
velocities are based on the velocities measured at the site and the values selected to 
represent the best estimate shear wave velocities are shown in Fig. 3-4. The best estimate 
input properties used in the equivalent linear analyses are listed in Appendix D. 


The modulus reduction curves shown in the upper part of Fig. 3-5 were used in this study. 
The upper curve was used for all the clay layers at the site; it is approximately equal to 


modulus reduction curve published by Sun et al (1988) and by Vucetic and Dobry (1991) 
for a plasticity index (PI) of about 30. The curve used for sand is approximately equal to 
the upper range for sand published by Seed and Idriss (1970). 


The damping ratio versus strain curve shown in the lower part of Fig. 3-5 was used for 
both the clay as well as the sand layers; it corresponds approximately to the lower range 
damping curve for sand published by Seed and Idriss (1970). 


3.3 Equivalent Linear Analyses 


The response of the Treasure Island site was calculated incorporating equivalent linear 
modulus and damping representation for the following conditions: , 


@ using the stronger component of the recording at Yerba Buena Island 
(YBI) as input rock outcrop motion, and 
¢ best estimate of the shear wave velocity profile 
¢ various estimates of the shear wave velocity profile 


® using the weaker component of the recording at Yerba Buena Island 
(YBI) as input rock outcrop motion, and 
¢ best estimate of the shear wave velocity profile 
¢ various estimates of the shear wave velocity profile 


@ using the best estimate of the shear wave velocity profile and the 
stronger and the weaker components of all the rock motions recorded in 
the San Francisco Bay Area (see Table B-3) as input outcrop motions 


The results of these analyses are presented below. 


3.3.1 Response at Treasure Island Using the Stronger Component of Recording 
at YBI as Input Rock Outcrop Motion 


3.3.1.1 Results Using Best Estimate Shear Wave Velocity Profile 


The best estimate shear wave velocities for the Treasure Island site are based on the 
measured shear wave velocities shown in Fig. 3-4 and are summarized below: 
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The corresponding values in SI units are listed below: 
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The response of this soil profile was calculated using the computer program SHAKE91. 
The results of the response calculations are presented in Figs. 3-6 through 3-9. Figure 3-6 
shows the calculated maximum shear strains and stresses and Fig. 3-7 shows the 
calculated strain-compatible damping ratios and shear wave velocities as well as the best 
estimate maximum shear wave velocities. As can be noted in these figures, the largest 
strain values are induced in the lower parts of the sandy fill and natural loose sand and in 
the Young Bay Mud layer. The maximum strains are quite small, however. Accordingly, 
the damping ratios are below 10% and the strain-compatible shear wave velocities are not 
much smaller than the maximum shear wave velocities as illustrated in Fig. 3-7. 


Figure 3-8 shows the variations of peak horizontal accelerations with depth. It is 
interesting that a significant part of the peak acceleration amplification occurs in the 
upper 40 feet (12.2 m) of the soil profile. 


The spectrum for the surface motion at Treasure Island, calculated using the best estimate 
shear wave velocities and the stronger component of the Yerba Buena Island record as 
input motion, is compared to the spectrum for the recorded surface motion in Fig. 3-9. 

As can be seen in this figure, the spectrum for the calculated surface motion provides a 


reasonable approximation of the spectrum for the recorded surface motion for almost the 
entire period range. 


Accordingly, the use of best estimate shear wave velocities and the recorded motion from 
a nearby rock site as input motion in an equivalent linear analysis provides a reasonable 
approximation for the recorded motions. 


The effects of using different shear wave velocities are considered in the next section. 
3.3.1.2 Results Using Various Shear Wave Velocity Profiles 


Several shear wave velocity profiles were considered for assessing the effects of varying 
shear wave velocities on the response calculations. These variations included using a 
lower bound and an upper bound shear wave velocity profiles and the addition of a 50-ft 
(15.2 m) layer of weathered shale at the bottom of the soil profile. They also included 
varying the shear wave velocity of the half-space from 3,000 to 5,000 ft/sec (914 to 1,524 
m/sec) and using two different values of the ratio of equivalent uniform strain to 
maximum strain. For most cases, a value of 0.5 was adopted for the this ratio, but for two 
of the cases a value of 0.35 was used; the latter value is similar to that used by Dickenson 
et al (1991) for similar calculations. 


The results of the response calculations, using various shear wave velocity profiles and 
the stronger component of the recording at Yerba Buena Island as input motion, are 
presented in Figs. 3-10, 3-11 and 3-12. Figure 3-10 shows the accelerograms calculated 
at the ground surface at the Treasure Island site for four of the cases considered. Also 
shown in the figure is the accelerogram of the stronger component of the recorded surface 
motion. As can be noted in the figure, the calculated motions are not significantly 
affected by the variations in the shear wave velocity profile. This is further illustrated in 
Fig. 3-11 which shows a comparison of the spectrum for the recorded motion and those 
for the calculated motions. Figure 3-12 indicates that the use of the best estimate shear 
wave velocities results in calculating a surface motion whose spectrum is very close to 
the median spectrum for all the cases considered. 


3.3.2 Response at Treasure Island Using the Weaker Component of the 
Recording at YBI as Input Rock Outcrop Motion 


3.3.2.1 Results Using Best Estimate Shear Wave Velocity Profile 
The same best estimate shear wave velocities which were used in conjunction with the 
stronger component of input motion were also used in conjunction with the weaker 


component of the input motion. 


The spectrum for the surface motion at Treasure Island, calculated using the best estimate 
shear wave velocities and the weaker component of the Yerba Buena Island record as 


input motion, is compared to the spectrum for the recorded surface motion in Fig. 3-13. 
The spectrum for the calculated surface motion is similar in shape to the spectrum for the 
recorded surface motion, but is considerably lower in amplitude over a significant range 
of periods. 


Accordingly, the use of best estimate shear wave velocities and the recorded motion from 
a nearby rock site as input motion in an equivalent linear analysis provides a reasonable 
approximation for the recorded motions only in one direction for the Treasure Island site. 


The effects of using different shear wave velocities are considered in the next section. 
3.3.2.2 Results using Various Shear Wave Velocity Profiles 


The same shear wave velocity profiles described in Section 3.3.1.2 were also considered 
for assessing the effects of varying shear wave velocities on the response calculations 
using the weaker component of the input motion. Figure 3-14 shows a comparison of the 
spectrum for the recorded motion and those for the calculated motions considering 
various shear wave velocity profiles and the weaker component of the Yerba Buena 
Island as input motion. Figure 3-15 shows the spectrum for the calculated surface motion 
using the best estimate shear wave velocities and the median spectrum for all the other 
cases considered together with the spectrum for the weaker component of the recorded 
surface motion. Use of the various shear wave velocity profiles appears to somewhat 
improve the comparison between the spectrum for the recorded motion and those for the 
calculated motions. 


3.3.3 Response at Treasure Island Using the Recorded Motions at All the Rock 
Sites in San Francisco Bay Area as Input Rock Outcrop Motions 


It is seldom that a soil site has an immediately adjacent rock site (such as was the case at 
Treasure Island and at Yerba Buena Island) both of which had instruments that recorded 
the motions generated by the same earthquake. Many times, recording stations at rock 
sites are several, if not tens, of kilometers away from the soil site under consideration. 


Several rock sites in the San Francisco Bay Area had strong motion instruments which 
were triggered during the Loma Prieta earthquake. Recordings at eleven such sites were 
obtained as summarized in Appendix B. These rock motions were then used as input 
rock outcrop motions for evaluating the response at the Treasure Island site. The best 
estimate shear wave velocities were used in all the response calculations. 


3.3.3.1 Results Using the Stronger Components 


The stronger component of each record obtained at the eleven rock sites was first scaled 
to a peak acceleration of 0.07g (1.e., its peak acceleration was made equal to that of the 
stronger component recorded at Yerba Buena Island). The response was then calculated 


using this scaled accelerogram as rock outcrop input motion using the best estimate shear 
wave velocities for the soil profile. Strain-compatible damping ratios and shear wave 
velocities were obtained for each response calculation. 


Figures 3-16a, 3-16b and 3-16c show a direct comparison of the spectrum for the 
recorded surface motion and the spectrum for the surface motion at Treasure Island 
calculated using the stronger component of motions recorded at each rock site. As can be 
seen from these figures, the results using the record at Yerba Buena Island provide the 
best estimate to the recorded values. The results using some of the other records (e.g., 
Piedmont and Sierra Point) provide reasonable estimates in some period ranges. The 
remaining records provide results that are significantly different from the recorded values. 


It may be noted that peak acceleration is estimated quite well using any of the rock 
records. Therefore, it is not sufficient to rely on making estimates of the peak 
acceleration in checking the accuracy of an analytical procedure. 


Figure 3-17 shows the spectrum for the recorded surface motion and the spectra for 
surface motions at Treasure Island calculated using the stronger component of motions 
recorded at Yerba Buena and at the other rock sites. This figure highlights the fact that 
using various rock input motions has a far greater effect on the calculated results than any 
other parameter. 


3.3.3.2 Results Using the Weaker Components 


The results using the weaker components are presented in Fig. 3-18. These results 
indicate a small improvement in the accuracy of estimating the spectrum of the recorded 
motion when the other rock records are used as input motions. 


3.3.3.3 Discussion of Results 


The material presented in this section of the report indicate that the results of a response 
analysis can be significantly affected by the input motion used in the response 
calculations. In fact, the input motion appears to have the most dominant effect on the 
results. Consequently, it would seem inappropriate to rely on the results of any one 
single analysis and that response calculations should always include variations in the 
input motion. 


In practice, it is quite often the case that the input motion is not known and may have to 
be generated using synthetically generated accelerograms or to rely on accelerograms 
recorded during prior earthquakes. The results shown in Figs. 3-16a, 3-16b, 3-16c, 3-17 
and 3-18 show the range of possible variations in computed spectral ordinates. Often, 
similar results are obtained for a particular site and then the median or the 84th percentile 
values of the ensemble of calculated spectra is used. This process is tested using the 


results for the Treasure Island site as summarized in Figs. 3-19 and 3-20. These figures 
show that: 


1. The median of the calculated spectra underestimates the spectrum for the recorded 
motion for either component or for the average of the two components. 


2. The 84th percentile spectrum for the calculated motions provides a reasonable 
estimate of the spectrum for the weaker component of the recorded motion, but does 
not do so for the stronger component as illustrated in Fig. 3-19. 


3. The 84th percentile spectrum for all the calculated motions (i.e., considering both 
components) appears to provide a reasonable estimate for the average spectrum for 
the two components of the recorded motions as shown in Fig. 3-20. 


3.4 Nonlinear Analyses 


The results using the program DESRA-?2 are presented in Fig. 3-21 together with those 
obtained using the results of an equivalent linear analysis. These results indicate that the 
DESRA-2? nonlinear analyses provide a somewhat reasonable estimate of the spectrum 
for the stronger component of motion at Treasure Island. They also show that the results 
of the equivalent linear analyses and those using this nonlinear procedure are comparable, 
although in this case the equivalent linear analysis provides an improved estimation of the 
spectrum of the recorded motion. Similar results were obtained using the nonlinear 
procedure incorporated in the program SUMDES. 


Similar comparisons have been made by Dickenson et al (1992) and by Chin and Aki 


(1993), who also concluded that at the levels of shaking experienced during the Loma 
Prieta earthquake, the two methods of analyses provide comparable results. 
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UNITED STATES OAKLAND WEST QUADRANGLE 
DEPARTMENT OF THE INTERIOR CALIFORNIA 
GEOLOGICAL SURVEY 7.5 MINUTE SERIES (TOPOGRAPHIC) 
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Fig. 3-1 Locations of Recording Stations at Treasure Island 
and at Yerba Buena Island 
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Fig. 3-2 Spectra for Surface Motions Recorded 
at Treasure Island (TI) and at Yerba Buena Island (YBI) 
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Fig. 3-3 Log of Boring and Shear Wave Velocities Measured at 
Treasure Island by USGS 
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Fig. 3-4 Shear Wave Velocities Measured at Treasure Island 
and Best Estimate Values of Velocities Used in the Analyses 
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Fig. 3-5 Modulus Reduction and Damping Relationships 
Used in All Analyses 
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Fig. 3-10 Accelerogram of Recorded Surface Motion and Accelerograms 
of Surface Motions at Tl Calculated Using Various Shear Wave Velocity 
Profiles and Stronger Component of YBI Record as Input Motion 
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4.0 SAN FRANCISCO AIRPORT SITE 


The recording station at San Francisco Airport is shown in Fig. 1-1 (Station No. 223) and 
in more detail in Fig. 4-1. The closest rock recording to the San Francisco Airport is that 
obtained at Sierra Point (Station No. 539 in Fig. 1-1). The recording station at Sierra 
Point is about 5 km north of the recording station at the San Francisco Airport. 


The spectra of the motions recorded at San Francisco Airport and at Sierra Point are 
presented in Fig. 4-2. The upper part of Fig. 4-2 shows the spectra for the stronger 
components of the recorded motions and the lower part of the figure shows the spectra for 
the weaker components. More details regarding the characteristics of the recordings at 
Sierra Point are given in Appendix B and those for San Francisco Airport are summarized 
in Appendix C. 


4.1 Subsurface Conditions at The San Francisco Airport Site 


The subsurface conditions at the San Francisco Airport site were obtained by drilling a 
boring nearby the station to a depth of 507 feet (154.5 m). The shear wave velocities 
were measured by USGS (Gibbs et al, 1992) and by Redpath (1991) using down-hole 
measuring techniques. 


The log of boring and the shear wave velocities measured by USGS are shown in Fig. 4-3 
and the shear wave velocities measured by USGS and by Redpath are shown in Fig. 4-4. 


The log of boring indicates that the subsurface conditions consist mostly of dense sands 
and Old Bay Mud to a depth of about 440 ft (134.1 m), except that there is a 15-ft (4.6 m) 
layer of Young Bay Mud at a depth of about 5 ft (1.5 m) below the ground surface. A 55- 
ft (16.8 m) layer of weathered mudstone exists below the soils and is underlain by a more 
component sandstone of the Franciscan formation as shown in Fig. 4-3. 


4.2 Dynamic Soil Properties 


The maximum shear wave velocities are based on the velocities measured at the site and 
the values selected to represent the best estimate shear wave velocities are shown in Fig. 
4-4. The best estimate input properties used in the equivalent linear analyses are listed in 
Subsection 4.3.1.1 below. 


The modulus reduction curves shown in the upper part of Fig. 3-4 and the damping ratio 


versus strain curve shown in the lower part of Fig. 3-4 were used in all the analyses 
conducted at the San Francisco Airport Site. 
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4.3 Equivalent Linear Analyses 


The response of the San Francisco Airport site was calculated incorporating equivalent 
linear modulus and damping representation for the following conditions: 


@ using the stronger component of the recording at Sierra Point (SP) as 
input rock outcrop motion, and 
¢ best estimate of the shear wave velocity profile 
¢ an upper and a lower range estimates of the shear wave 
velocity profile 


@ using the weaker component of the recording at Sierra Point (SP) as 
input rock outcrop motion, and 
¢ best estimate of the shear wave velocity profile 
¢ an upper and a lower range estimates of the shear wave 
velocity profile 


@ using the best estimate of the shear wave velocity profile and the 
stronger and the weaker components of all the rock motions recorded in 
the San Francisco Bay Area (see Table B-3) as input outcrop motions 


The results of these analyses are presented below. 


4.3.1 Response at San Francisco Airport Using the Stronger Component of 
Recording at Sierra Point as Input Rock Outcrop Motion 


4.3.1.1 Results Using Best Estimate Shear Wave Velocity Profile 


The best estimate shear wave velocities for the San Francisco Airport site are based on 
the measured shear wave velocities shown in Fig. 4-4 and are summarized below: 
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The corresponding values in SI units are listed below: 
Avproximate Depth Below 
Ground Surface - meter 

384.1 


Maximum Shear 
Wave Velocity 
m/sec 


Total Unit 
Weight 
kN/m* 


samara 3 


The response of this soil profile was calculated using the computer program SHAKE91. 
The results of the response calculations are presented in Figs. 4-5 through 4-8. Figure 4-5 
shows the calculated maximum shear strains and stresses and Fig. 4-6 shows the 
calculated strain-compatible damping ratios and shear wave velocities. As can be noted 
in these figures, the largest strain values are induced in the Young Bay Mud layer. The 
maximum strains are quite small, however. Accordingly, the damping ratios are not very 
high and except for the Young Bay Mud layer the strain-compatible shear wave velocities 
are only slightly lower than the maximum shear wave velocities listed above. 


Figure 4-7 shows the variations of peak horizontal accelerations with depth. It is 
interesting that a significant part of the peak acceleration amplification occurs in the 
upper 20 ft (6.1 m) of the soil profile. 


The spectrum for the surface motion at San Francisco Airport, calculated using the best 
estimate shear wave velocities and the stronger component of the Sierra Point record as 
input motion, is compared to the spectrum for the recorded surface motion in Fig. 4-8. 
The spectrum for the calculated surface motion provides a very good approximation of 
the spectrum for the recorded surface motion for almost the entire period range. 


The effects of using different shear wave velocities are considered in the next section. 
4.3.1.2 Results using Various Shear Wave Velocity Profiles 
In addition to the best estimate shear wave velocities, an upper range and a lower range 


shear wave velocity profiles were considered for assessing the effects of varying shear 
wave velocities on the response calculations at the San Francisco Airport site. 


The results of the response calculations, using these three shear wave velocity profiles 
and the stronger component of the recording at Sierra Point as input motion, are presented 
in Fig. 4-9. Figure 4-9 shows a comparison of the spectrum for the recorded motion and 
those for the calculated motions. These results indicate that the use of the lower range 
values of shear wave velocities appear to underestimate the spectral ordinates for the 
recorded motion for periods shorter than about 0.4 sec. The use of the upper range values 
of shear wave velocities appears to overestimate these spectral ordinates in the period 
range of about 0.3 to 0.5 sec. The use of all three sets of shear wave velocities, however, 
appears to underestimate the peak acceleration. 


4.3.2 Response at San Francisco Airport Using the Weaker Component of the 
Recording at Sierra Point as Input Rock Outcrop Motion 


The same three sets of shear wave velocities (i.e., lower range, best estimate and upper 
range shear wave velocities) that were used in conjunction with the stronger component 
of input motion were also used in conjunction with the weaker component of the input 
motion. 


The results of the response calculations are presented in Fig. 4-10. The spectra for the 
calculated surface motions are similar in shape to the spectrum for the recorded surface 
motion, but are considerably lower in amplitude over a significant range of periods. 


Accordingly, the use of best estimate shear wave velocities and the recorded motion from 
a nearby rock site as input motion in an equivalent linear analysis appears to provide a 
reasonable approximation for the recorded motions only in one direction for the San 
Francisco Airport site. This finding is similar to that obtained for the Treasure Island site 
(Section 3.3.2.1) 


4.3.3 Response at San Francisco Airport Using the Recorded Motions at All the 
Rock Sites in San Francisco Bay Area as Input Rock Outcrop Motions 


The stronger component of each record obtained at the eleven rock sites was first scaled 
to have a peak acceleration of 0.1g (1.e., its peak acceleration was made equal to that of 
the stronger component recorded at Sierra Point). The response was then calculated using 
this scaled accelerogram as rock outcrop input motion using the best estimate shear wave 
velocities for the soil profile. Strain-compatible damping ratios and shear wave velocities 
were obtained for each response calculation. 


Figure 4-11 shows the spectrum for the recorded surface motion and the spectra for 
surface motions at San Francisco Airport calculated using the stronger component of 
motions recorded at Sierra Point and at the other rock sites. This figure highlights the 
fact that using various rock input motions has a far greater effect on the calculated results 
than any other parameter. 
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Similar calculations were completed using the weaker components of the records 
obtained at all the rock sites in the San Francisco Bay Area. The weaker component of 
each record obtained at the eleven rock sites was first scaled to have a peak acceleration 
of 0.06g (i.e., its peak acceleration was made equal to that of the weaker component 
recorded at Sierra Point). 


The results using the weaker components are presented in Fig. 4-12. These results 
indicate an improvement in the accuracy of estimating the spectrum of the recorded 
motion when the other rock records are used as input motions. 


As was the case for the Treasure Island site, the results of a response analysis can be 
significantly affected by the input motion used in the response calculations. In fact, the 
input motion appears to have the most dominant effect on the results. Consequently, it 
would seem inappropriate to rely on the results of any one single analysis and that 
response calculations should always include variations in the input motion. 


Again the use of the median or the 84th percentile values of the ensemble of calculated 
spectra is tested using the results for the San Francisco Airport site as summarized in 
Figs. 4-13 and 4-14. These figures show that: 


1. The median of the calculated spectra underestimates the spectrum for the recorded 
motion for either component or for the average of the two components. 


2. The 84th percentile spectrum for the calculated motions provides a reasonable 
estimate of the spectrum for both components of the recorded motion, except at a 
narrow period range as illustrated in Fig. 4-13. 


3. The 84th percentile spectrum for all the calculated motions (i.e., considering both 
components) appears to provide a reasonable estimate for the average spectrum for 
the two components of the recorded motions except over a very narrow period range 
as shown in Fig. 4-14. 
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Fig. 4-1 Location of Recording Station at The 
San Francisco Airport 
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Fig. 4-2 Spectra for Surface Motions Recorded 
at the San Francisco Airport (SFO) and at Sierra Point (SP) 
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Fig. 4-3 Log of Boring and Shear Wave Velocities Measured at 
the San Francisco Airport by USGS 
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Fig. 4-4 Shear Wave Velocities Measured at The San Francisco 
Airport and Best Estimate Values of Velocities Used in the Analyses 
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5.0 SUMMARY AND CONCLUSIONS 


The Loma Prieta earthquake was the largest earthquake to affect the San Francisco Bay 
Area since the 1906 San Francisco earthquake. It triggered by far the largest number of 
strong motion recording stations ever triggered by an earthquake. It provided earthquake 
ground motions at many free field stations underlain by various subsurface conditions. 


The peak accelerations recorded at all these free field stations are summarized in 
Appendix A. The characteristics of the horizontal earthquake ground motions recorded at 
rock sites are summarized in Appendix B and those related to the horizontal motions 
recorded at soft soil sites are summarized in Appendix C. 


The characteristics summarized in these appendices provide a significant insight into the 
nature of earthquake ground motions and the potential effects of site conditions and 
distance on these characteristics. 


For example, the data indicate that peak horizontal accelerations at rock sites are 
comparable to those recorded at soil sites (other than soft soil sites). The peak horizontal 
accelerations at soft soil sites, however, are significantly larger than those at the other two 
site conditions. 


The spectral shapes for horizontal motions recorded at rock sites show a strong distance 
dependence as shown in Fig. B-39. Typically, most available attenuation relationships 
assume a distance-independent shape or a shape that is weakly dependent on distance. 
This trend is illustrated in Fig. 5-1 which provides a comparison of the median spectral 
shapes for horizontal rock motions recorded during the Loma Prieta earthquake (Fig. B- 
39) with spectral shapes calculated using the attenuation relationship for rock sites 
derived by Idriss (1991b). This attenuation relationship incorporates distance- 
dependence of spectral shapes, but apparently not sufficiently to account for the measured 
effects obtained from Loma Prieta. Incorporation of the distance dependence of spectral 
shapes is essential for developing appropriate target rock spectra for a specific application 
and for conducting ground response analyses. 


The records obtained at the soft soil sites offered an excellent opportunity to assess the 
procedures used for calculating response of soil sites during earthquakes. This was 
particularly the case because rock outcrop motions were also available from nearby 
locations and subsurface information was gathered at most of these soft soil sites. 


Examination and studies of the response at these soft soil sites have been completed by 
several investigators since the occurrence of the Loma Prieta earthquake. These 
investigators include Borcherdt and Glassmoyer (1991), Chin and Aki (1993), Dickenson 
et al (1991), EPRI (1993), Hryciw et al (1991), Idriss (1990, 1991a), Jarpe et al (1989), 
Rollins et al (1992), Schneider et al (1991) and Yokel (1992). These investigations 
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included examination of the recorded values at the soft soil sites in comparison to those 
recorded at the rock sites. 


The results of the studies included in this report provide the following observations: 


I 


The equivalent linear procedure for calculating response of soil sites can provide 
reasonable estimates of the recorded motions if dynamic soil properties are 
reasonably well defined and if the rock outcrop input motion is reasonably well 
defined. 


The estimation of earthquake ground motions using ground response analyses is 
affected greatly by the input motion. Therefore, it is important that variations in this 
input motion be incorporated in estimating earthquake ground motions for future 
events. 


It is recommended that the 84th percentile of the results of site response analyses be 
used as an estimate of future shaking levels for the earthquake magnitude and 
distance under consideration. 


It is important that response calculations take into account the difference in frequency 
content of rock motions recorded at close distances to the source compared to the 
frequency content of rock motions recorded at larger distances. This is particularly 
important when conducting site response analyses for sites located close to the 
potential source of future earthquakes. 


Future studies need to address the response of deep soil sites (other than soft soil sites) 
and to make a systematic examination of the factors that affect the response of such sites. 
Subsurface data at many of these sites have been recently collected and can be used in 
studies similar to those presented in this report. 
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APPENDIX A 


EARTHQUAKE GROUND MOTIONS RECORDED 
DURING THE 1989 LOMA PRIETA EARTHQUAKE 


To Report on: 


ASSESSMENT OF SITE RESPONSE 
ANALYSIS PROCEDURES 


APPENDIX A 


EARTHQUAKE GROUND MOTIONS RECORDED 
DURING THE 1989 LOMA PRIETA EARTHQUAKE 


A.1 INTRODUCTION 


The Loma Prieta earthquake occurred on October 17, 1989 at 5:04 pm Pacific daylight 
time along a 45-km long segment of the San Andreas fault in the Santa Cruz Mountains 
as shown in Fig. A-1. The earthquake was named after the highest topographic point 
(3791 ft or approximately 1155 m) ) adjacent to the fault zone. The hypocenter of the 
earthquake was at a depth of about 18 km. The rupture plane dips to the southwest at 
about 70 degrees, as shown in Fig. A-2; thus the epicenter is several kilometers west of 
the San Andreas fault trace. The earthquake was assigned a surface wave magnitude, M, 
= 7.1, and a moment magnitude, M,, = 7. The rupture was bilateral, i.e., the rupture was 
initiated at the hypocenter and propagated both northward and southward simultaneously. 


Strong motion instruments had been installed by the United States Geological Survey 
(USGS) and by the California Strong Motion Instrumentation Program (CSMIP) at 
numerous locations and on various site conditions. These included: 33 stations at rock 
sites (of these two stations were at abutments of earth dams); nine stations at soft soil 
sites; and 47 stations at other soil sites (other than soft soil sites). The locations of the 
USGS stations that recorded the Loma Prieta earthquake are shown in Fig. A-3 and the 
locations of the CSMIP stations are shown in Fig. A-4. 


The peak horizontal accelerations and the peak vertical accelerations recorded during this 
earthquake at free-field sites are summarized below. 


The characteristics of the horizontal motions recorded at rock sites are summarized in 
Appendix B and the characteristics of the horizontal motions recorded at the soft soil sites 
are summarized in Appendix C. 


A.2 PEAK ACCELERATIONS OF RECORDED MOTIONS 


The Volume I" peak accelerations (two horizontal and one vertical components) of the 
motions recorded by the USGS and by CSMIP at rock sites are listed in Table A-1. The 
Volume I peak accelerations of the motions recorded at soft soil sites are listed in Table 


Volume I values are digitized at an unequal time basis, aimed at defining all significant 
features of the record (Hudson, 1976), and hence contain all the peak values of the original 
record. Volume II values, however, are digitized at an equal time increment and can result in 
peak accelerations that are somewhat smaller than the peak values in the original record. 
Volume II records (i.e., those digitized at an equal time increment) are usually used to calculate 
velocities, displacements, response spectra, Fourier amplitudes ... etc. 
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A-2 and those recorded at other soil sites are listed in Table A-3. The data listed in these 
tables are for stations considered to represent free-field conditions (except possibly for 
those at the abutments of earth dams and at the Emeryville site). 


Included in each table are the name of the station, the distance to the source and the 
volume I peak horizontal and vertical accelerations. Note that the distance given in these 
tables pertains to the distance from the recording station to the closest point along the 
rupture surface (as defined by aftershocks) at a depth of three kilometers below the 
ground surface. 


The horizontal peak accelerations are shown in Fig. A-5. The upper part of the figure 
shows the peak horizontal accelerations of both the strong and the weak components and 
the lower part of the figure shows the geometric average of the two horizontal 
components. As can be noted in both parts of the figure, the peak horizontal accelerations 
at rock and at soil sites (other than soft soil sites) are very similar and little distinction can 
be made between the two sets of peak horizontal accelerations. The peak horizontal 
accelerations at the soft soil sites, however, are significantly greater than those at the 
other two site conditions. 


The vertical accelerations are shown in the upper part of Fig. A-6. The lower part of the 
figure shows the ratio of the peak vertical acceleration, a,, divided by the geometric 
average, a,, of the two horizontal components. The peak vertical accelerations at rock 
sites appear to be smaller than those recorded at soil sites during this earthquake. The 
variations of the ratio a,/ a; with distance for the data recorded during this earthquake are 
typical of such variations. The ratio for soft soil sites is quite low (because of the high 
horizontal accelerations) and the data points in Fig. A-6 are among the few available for 
this ratio at soft soil sites. 
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Table A-1 
Volume I Peak Accelerations of Earthquake Ground Motions Recorded by USGS and by 
CSMIP at Rock Sites During the Loma Prieta Earthquake 


0.640 | 0500 | 0.470 
Gilroy #1 10.0 0.500 0.430 0.220 


Gilroy #6 0.170 0.130 0.100 


Abutment - Coyote Lake Dam 0.490 0.150 0.080 
0.470 0.440__[ 0.400 
Hollister (S St. & Pine Drive) 0.060 0.040 0.050 


Sago South (Cienega Road) 0.0705 art ume Uae aa 
Woodside Fire Station 0.080 0.080 0.050 


Abutment - Lexington Dam 0.450 0.410 0.150 
Redwood City 0.090 0.050 0.040 


Monterey - City Hall 0.070 0.070 0.030 


0.070 0.080 


0,060 
0,060 
0.040 
52.0 0,050 
0.030 
0,050 
0.030 
0.050 
0.030 
0.030 
0.030 
0.030 
508 0.030 
0.030 
0.020 
0.060 
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Table A-2 
Volume I Peak Accelerations of Earthquake Ground Motions Recorded by USGS and by 
CSMIP at Soft Soil Sites During the Loma Prieta Earthquake 


[Larkspur Ferry Terminal | 99.0| 0.140 | 0.100 | 0.060 


¢¢ Peak accelerations for this site are Volume II values. 
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Table A-3 
Volume I Peak Accelerations of Earthquake Ground Motions Recorded by USGS and by 
CSMIP at Soil Sites (other than Soft Soil Sites) During the Loma Prieta Earthquake 


0.270 
0.380 
0.170 
0.600 


Anderson Dam - downstream 0.260 0.250 0.170 
Coyote Lake Dam - downstream 0.190 0.170 0.100 
San Jose (Santa Teresa Hills) 0.280 0.270 0.220 


Gilroy #7 0.330 0.100 0.120 
Hollister Airport 0.290 0.270 0.160 


0.190 0.100 


0.100 
[Sunnyvale - Colton Ave. | 27.1_| 0.220 | _0.190_| 

0.220 
0.200 


0.080 
0.200 
0,060 
0.110 
0.090 
0.070 
0.100 
0.110 
0,090 
0.070 
0.030 
0.100 
0.100 
0.080 
0.040 
0.060 
0.030 
0.020 


San Jose (Station No. 57562) 0.200 0.200 0.140 
Halls Valley - Grant Park 0.130 0.060 0.110 
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Table A-3 (Cont'd) 
Volume I Peak Accelerations of Earthquake Ground Motions Recorded by USGS and by 
CSMIP at Soil Sites (other than Soft Soil Sites) During the Loma Prieta Earthquake 
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BB = Bay Bridge 
NF = Nimitz Freeway 
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Fig. A-1 Location of the Main Shock of the Loma Prieta Earthquake, Area 
of Aftershocks and Inferred Extent of Fault Rupture (from McNutt, 1990) 
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Fig. A-2. Transverse Section Showing 
Distribution of Aftershocks and Dip of 
Fault Plane (from McNutt, 1990) 
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Fig. A-3 Locations of USGS Stations Which 


Recorded the Loma Prieta Earthquake 
(from Maley et al, 1989) 
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Fig. A-4 Locations of CSMIP Stations Which 


Recorded the Loma Prieta Earthquake 
(from Shakal et al, 1989) 
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Fig. A-5 Horizontal Accelerations Recorded During the 1989 Loma Prieta Earthquake 
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Fig. A-6 Vertical Accelerations Recorded During the 1989 Loma Prieta Earthquake 
and Ratio of Vertical Acceleration to Horizontal Acceleration 
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CHARACTERISTICS OF THE HORIZONTAL EARTHQUAKE GROUND 
. MOTIONS RECORDED AT ROCK SITES DURING 
THE 1989 LOMA PRIETA EARTHQUAKE 


B.1 INTRODUCTION 


As described in Appendix A, 33 of the 89 free field recordings obtained during Loma 
Prieta earthquake were at rock sites. Of particular interest to this study are the recordings 
obtained in the vicinity of the source (say, within 20 + km) and those obtained in the San 
Francisco-Bay Area. The peak accelerations, peak velocities, peak displacements, 
spectral ordinates and effective duration (as defined from a Husid plot) of these rock 
motions are summarized in this Appendix. 


B.2 ROCK MOTIONS WITHIN 20 KM FROM THE SOURCE 


Six recording stations located on rock were within 20 km from the source of the Loma 
Prieta earthquake. Two of these recording stations were at abutments of earth dams; the 
characteristics of the ground motions obtained from these two stations are presented in 
this Appendix but are not used in evaluating the characteristics of rock motions within 20 
km of the source. 


The characteristics of the horizontal earthquake ground motions recorded at these rock 
sites during the Loma Prieta earthquake are listed in Tables B-1 through B-3 and in Figs. 
B-1 through B-14 and are summarized below. 


B.2.1 Peak Horizontal Accelerations, a, Velocities, v, and Displacements, d 


The peak accelerations (Volume II), velocities and displacements of the earthquake 
ground motions recorded at these rock sites are listed in Table B-1. 


B.2.2 Ratios v/a and ad/v2 


The ratios v/a and ad/v? for the twelve horizontal components listed in Table B-1 are 
listed in Table B-2. The average value of the ratio v/a for all the records listed in Table 
B-1 is about 110 cm/sec/g if the records obtained at the dam abutments are included. If 
these two records are not included the average value of the ratio v/a is about 79. The 
corresponding average values of the ratio ad/v2 are about 2.9 and 3.5, respectively. 
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B.2.3 Accelerograms and Spectral Accelerations 


The accelerograms and the spectral accelerations (5 percent spectral damping) for these 
ground motions are presented in Figs. B-1 through B-6. The normalized spectra for the 
horizontal accelerograms recorded at Corralitos, Gilroy No. 1, Gilroy No. 6 and Santa 
Cruz (i.e., normalized spectra for all components listed in Table B-1 except for the four 
components recorded at abutments of earth dams) are shown in Fig. B-7. Also shown in 
Fig. B-7 is the median normalized spectrum for these eight accelerograms. 


B.2.4 Husid Plots and Effective Duration 


Effective duration of a given accelerogram is assessed using the procedure originally 
proposed by Trifunac and Brady (1975). This procedure is based on the time required for 
the buildup of the integral 


ee (t dt 


0 


in which a(t) is the acceleration time history. Arias (1969) showed that this integral is a 
measure of the energy of the accelerogram, and defined the intensity of the entire record 
by the following expression: 


tf 
= 5g Jar (that 


in which J, is the Arias’ intensity and tfis the total duration. Husid (1969) proposed the 
use of the normalized variable h(t): 


L(t) la°(1 Jat 


i) 


tf Equation B-1 
\a2(t at 

0 
Thus, h(t) = 0 at the beginning of the record and = 1 (or 100 percent) at the end of the 


record. The plot of h(t) versus ¢ is designated the Husid plot. 


Effective duration as defined by Trifunac and Brady (1975) is the time interval needed for 
h(t) to build up from 5 to 95 percent. 


The Husid plots for the twelve horizontal components listed in Table B-1 are presented in 
Figs. B-8 through B-13. The Husid plots for all the motions considered are shown in Fig. 
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B-14. The time at which h(t) = 5 percent and the time at which h(t) = 95 percent are 
listed in Table B-3 for each component included in Figs. B-8 through B-13. Also listed 
in Table B-3 are the values of effective duration as defined above. The effective duration 
varies from about 3.7 to 8 seconds within distances of about 10 km from the source and 
from about 9.5 to 13 seconds at distances of 18 to 20 km (the range is from 9.5 to about 
16 seconds if the recordings at the abutment of the Coyote dam are included). 


B.3 ROCK MOTIONS IN THE SAN FRANCISCO BAY AREA 


The recording stations located on rock sites in the San Francisco Bay Area are listed in 
Table B-4. The characteristics of the horizontal earthquake ground motions recorded at 
these eleven rock sites during the Loma Prieta earthquake are listed in Tables B-4 through 
B-6 and in Figs. B-15 through B-38 and are summarized below. 


B.3.1 Peak Horizontal Accelerations, a, Velocities, v, and Displacements, d 


The peak accelerations (Volume II), velocities and displacements of the earthquake 
ground motions recorded at these rock sites are listed in Table B-4. 


B.3.2 Ratios v/a and ad/v2 


The ratios v/a and ad//? for the twelve records listed in Table B-4 are listed in Table B- 
5. The average value of the ratio v/a for all the records listed in Table B-4 is about 147 
cm/sec/g and the average value of ad/v? is about 2.1. 


B.3.3 Accelerograms and Spectral Accelerations 


The accelerograms and the spectral accelerations (5 percent spectral damping) for these 
ground motions are presented in Figs. B-15 through B-25. The normalized spectra for all 
these horizontal accelerograms are shown in Fig. B-26. Also shown in Fig. B-26 is the 
median normalized spectrum for these 22 accelerograms. 


B.3.4 Husid Plots and Effective Duration 


The Husid plots (using Eq. B-1) for the horizontal accelerograms obtained at these rock 
sites in the San Francisco Bay Area are presented in Figs. B-27 through B-37. The Husid 
plots for the motions considered are shown in Fig. B-38; note that the Husid plots for the 
two horizontal components at Telegraph Hill are not included in Fig. B-38 because 
several seconds of the early part of these two components had been eliminated in the 
digitizing process. The time at which h(t) = 5 percent and the time at which h(t) = 95 
percent are listed in Table B-6 for each component included in Figs. B-27 through B-37. 
Also listed in Table B-6 are the values of effective duration as defined above; the values 
of effective duration range from about 6 to 20 seconds. 
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B.4 DISCUSSION 


The material presented in this Appendix provides information regarding the general 
characteristics of rock motions recorded at close distance from the source and those 
recorded at considerably larger distances from the source. The median frequency content 
of each set of recordings is shown in Fig. B-39 in terms of normalized spectral shapes. 

As can be noted in the figure, the median normalized spectrum for the records obtained at 
close distances from the source has a peak at a period of about 0.22 sec (or a frequency of 
about 4.5 Hz) while the median normalized spectrum for the more distant records has a 
peak at about a period of about 0.42 sec (or a frequency of about 2.4 Hz). These 
frequency characteristics are very important to take into account when selecting input 
motions for response calculations, especially for softer soil sites (Idriss, 1991). 


The effective duration values obtained for all the recordings on rock presented in this 
Appendix are shown in Fig. B-40. Also shown in the figure is the line representing the 
least square fit for the data shown in the figure; the equation of this line is as follows: 


Ln(D,) = 1.695 + 0.151Ln(R) Equation B-2 


in which D., is the effective duration in seconds and R is distance (as defined in 
Appendix A) in km. The standard error term is 0.31 (natural logarithm basis). The trend 
represented by the data and the results of the least square fit suggest a weak dependence 
of effective duration on distance. 
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Table B-1 
Peak Accelerations, Velocities and Displacements of Horizontal Earthquake Ground 
Motions Recorded at Rock Sites Within a Distance of about 20 km from the Source 
during the Loma Prieta Earthquake 


Peak Peak 
seg Acceleration | Velocity 
Station Component g cm/sec 


ee 0478 
[NS [0.629 a 


Peak 
ae 


sep 
TNS TP 0.44a gad 
Giioy Fi Re irh 2 a ea 442 
PEUNS 2" [ 0435 ;roy = 319" n |e al 

Gilroy #6 18.1 


Abutment - Coyote | 19.5 
he bat ON ae 
Notes: 

1. Peak accelerations listed in this table are "Volume II" accelerations. 
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Table B-2 
Values of v/a and ad / v’ of Horizontal Earthquake Ground Motions Recorded at 
Rock Sites Within a Distance of about 20 km from the Source during the Loma Prieta 
Earthquake 


Distance via 
Spe 
MMM 22 2a 
Abutment - Lexington SSEWs—t 
Shae RE EE Se ere aes 
aaa == Sa! 


Abutment - Coyote [NSE 
Lake Dam =e 
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Table B-3 
Effective Duration Determined from Husid Plots of Horizontal Earthquake Ground Motions 
Recorded at Rock Sites Within a Distance of about 20 km from the Source during the Loma 


Prieta Earthquake 
ere Time at 5% | Time at 95% Duration 
Station Component seconds seconds seconds 
De saa iewe 
os na ee 
Abutment ee 
Gilroy #1 


Effective 


Gilroy #6 18.1 
a Eo 

19.5 
Lake Dam 


Santa Cruz 20.2 13.95 
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Table B-4 
Peak Accelerations, Velocities and Displacements of Horizontal Earthquake Ground Motions 
Recorded at Rock Sites in San Francisco-Bay Area during the Loma Prieta Earthquake 


me 
Station Component 


Sierra Point 
ee 
Diamond Heights | 76 
a 
Rincon Hill 75 

ie SSE eee 
Pacific Heights 
8 Se Se 
Telegraph Hi 81 
le 
Golden Gate 82 
22) SS 

Presidio 82 
Pie oe? 
Cliff House 33 

2 SCS er ee ee ee 

Yerba Buena 79 

Se Se eee 


Peak 
cereale i 


Acceleration |Peak Velocity 


Piedmont Nasw | __0.083 
N&SE | 0.071 


Lawrence Berkeley | 83 0.117 
Laboratory 0.049 


Note: peak accelerations listed in this table are "Volume II" accelerations. 
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Table B-5 
Values of v/ a and ad / v’ of Horizontal Earthquake Ground Motions Recorded at Rock Sites 
in San Francisco-Bay Area during the Loma Prieta Earthquake 


heal v/ a 
Station Component | cm/sec/g zh 
Sierra Point se 
fae 
Diamond Heights Fre eae cael 
ee soe 
Rincon Fil PEW. aloe sooe 
Bea Coed rte eer es 


Pacific Heights See ee cee eee 
[5 NS a eipe1Oh | eld oba 
Telegraph il | 81 


82 
Presidio 82 
oF RS NS 3 anole 
Cliff House 
ee ao 
Yerba Buena 
com aman 
Piedmont eS ae 


Lawrence Berkeley 
Laboratory 


Page B-9 


Table B-6 
Effective Duration Determined from Husid Plots of Horizontal Earthquake Ground Motions 
Recorded at Rock Sites in San Francisco-Bay Area during the Loma Prieta Earthquake 


Effective 

Distance Time at 5% | Time at95%]| Duration 

seconds seconds seconds 
Sierra Point 
ee arts 
Diamond Heights | 76 
eg eee ete 
Rincon Hill 79 
ee eo 
Pacific Heights 
al 
Telegraph Hil a1 
ene ne 
Golden Gate [ew [9.76 [iss | 6.05 
ee ee es ee 
Presidio AEP YT aS PCT 
SS 
Cliff House 33 
ea 
Yerba Buena 79 
OR eee ale 


Piedmont N45W 18.96 11.76 
N45E 19.11 11.70 


Lawrence Berkeley MeEW 17.14 
ce a Sie eee waaay Ne 
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Fig. B-1 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Corralitos during the Loma Prieta Earthquake 
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Fig. B-2 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Abutment of Lexington Dam during the Loma Prieta Earthquake 
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Fig. B-3 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Gilroy No.1 during the Loma Prieta Earthquake 
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Fig. B-4 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Gilroy No. 6 during the Loma Prieta Earthquake 
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Fig. B-5 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Abutment of Coyote Dam during the Loma Prieta Earthquake 
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Fig. B-6 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Santa Cruz during the Loma Prieta Earthquake 
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Fig. B-7 Normalized Spectra for Motions Recorded at Rock 
Sites within 20 km from Source During the Loma Prieta Earthquake 
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Fig. B-8 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Corralitos during the Loma Prieta Earthquake 
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Fig. B-9 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Abutment of Lexington Damduring the Loma Prieta Earthquake 
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Fig. B-10 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Gilroy No.1 during the Loma Prieta Earthquake 
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Fig. B-11 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Gilroy No. 6 during the Loma Prieta Earthquake 
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Fig. B-12 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Abutment of Coyote Dam during the Loma Prieta Earthquake 
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Fig. B-13 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Santa Cruz during the Loma Prieta Earthquake 
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Fig. B-15 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Sierra Point during the Loma Prieta Earthquake 
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Fig. B-16 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Diamond Heights during the Loma Prieta Earthquake 
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Fig. B-17 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Rincon Hill during the Loma Prieta Earthquake 
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Fig. B-18 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Pacific Heights during the Loma Prieta Earthquake 
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Fig. B-19 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Telegraph Hill during the Loma Prieta Earthquake 
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Fig. B-20 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Golden Gate during the Loma Prieta Earthquake 
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Fig. B-21 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at the Presidio during the Loma Prieta Earthquake 
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Fig. B-22 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at the Cliff House during the Loma Prieta Earthquake 
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Fig. B-23 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Yerba Buena Island during the Loma Prieta Earthquake 
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Fig. B-24 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Piedmontd during the Loma Prieta Earthquake 
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Fig. B-25 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Lawrence Berkeley Lab during the Loma Prieta Earthquake 
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Fig. B-26 Normalized Spectra for Motions Recorded at Rock Sites 
in the San Francisco Bay Area during the Loma Prieta Earthquake 
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Fig. B-27 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Sierra Point during the Loma Prieta Earthquake 
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Fig. B-28 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Diamond Heights during the Loma Prieta Earthquake 
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Fig. B-29 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Rincon Hill during the Loma Prieta Earthquake 
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Fig. B-30 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Pacific Heights during the Loma Prieta Earthquake 
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Fig. B-31 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Telegraph Hill during the Loma Prieta Earthquake 
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Fig. B-32 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Golden Gate during the Loma Prieta Earthquake 
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Fig. B-33 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at the Presidio during the Loma Prieta Earthquake 
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Fig. B-34 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at the Cliff House during the Loma Prieta Earthquake 
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Fig. B-35 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Yerba Buena Island during the Loma Prieta Earthquake 
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Fig. B-36 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Piedmont during the Loma Prieta Earthquake 
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Fig. B-37 Husid Plots of Horizontal Components of Earthquake Ground Motions 
Recorded at Lawrence Berkeley Laboratory during the Loma Prieta Earthquake 


eyenbypylez ejalid eulo7] ay} Buiing eaiy Aeg oosiouesy UeS ay} UI SazIS YIOY }e 
paepsorsay SuoHop Ppunosy ayenbyplez jo syuauodwioy jeJUuOZUOH Jo sjojd pisny g¢E-g ‘big 


2aS - sul] 


D fi iid 02 


| il 


OV 


Ni 
CP } 09 


iy 
My 


08 


OOF 


[We,QJe] pazijeuuoy eane;nuing 


Spectral Acceleration / Zero Period Acceleration 


Median Spectral 

V4 Shape for Rock 
YA—| Motions Recorded 
in San Francisco 
Bay Area (Fig. B-26) 


“0.01 0.10 1.00 


Period - sec 


Fig. B-39 Median Normalized Spectral Shapes for Rock 
Motions Recorded During the Loma Prieta Earthquake 
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Fig. B-40 Effective Duration Obtained for Motions 
Recorded at Rock Sites During the Loma Prieta Earthquake 
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APPENDIX C 


CHARACTERISTICS OF THE HORIZONTAL 
EARTHQUAKE GROUND MOTIONS RECORDED AT 
SOFT SOIL SITES DURING THE 1989 LOMA PRIETA 

EARTHQUAKE 


To Report on: 


ASSESSMENT OF SITE RESPONSE 
ANALYSIS PROCEDURES 


APPENDIX C 


CHARACTERISTICS OF THE HORIZONTAL EARTHQUAKE GROUND 
MOTIONS RECORDED AT ROCK SITES DURING THE 1989 LOMA 
PRIETA EARTHQUAKE 


C.1 INTRODUCTION 


As described in Appendix A, nine of the 92 free field recordings obtained during Loma 
Prieta were at soft soil sites, all of which are in the San Francisco-Bay Area (see Fig. 1). 
The peak accelerations, peak velocities and peak displacements, spectral ordinates and 
effective duration (as defined from a Husid plot) of these rock motions are summarized in 
this Appendix. 


C.2 CHARACTERISTICS OF HORIZONTAL MOTIONS RECORDED AT SOFT 
SOIL SITES 


The characteristics of the horizontal earthquake ground motions recorded at these nine 
soft soil sites during the Loma Prieta earthquake are listed in Tables C-1 through C-3 and 
in Figs. C-1 through C-21 and are summarized below. 


C.2.1 Peak Horizontal Accelerations, a, Velocities, v, and Displacements, d 


The peak accelerations (Volume II), velocities and displacements of the earthquake 
ground motions recorded at these rock sites are listed in Table C-1. 


C.2.2 Ratios v/a and ad/v2 


The ratios v/a and ad/y for the twelve records listed in Table C-1 are listed in Table C- 
2. The average value of the ratio v/a for all the records listed in Table C-1 is about 147 
cm/sec/g and the average value of the ratio ad/v2 is about 1.6. Note that the recordings at 
Emeryville were not used to obtain the average values of v/a and ad/v?; the effect of 
including the values for these recordings, however, is minimal. 


C.2.3 Accelerograms and Spectral Accelerations 
The accelerograms and the spectral accelerations (5 percent spectral damping) for these 


ground motions are presented in Figs. C-1 through C-9. The normalized spectra for all 
the recordings, except for those at Emeryville, are shown in Fig. C-10. 
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C.2.4 Husid Plots and Effective Duration 


The Husid plots (using Eq. B-1 in Appendix B) for the horizontal accelerograms obtained 
at these soft soil sites in the San Francisco Bay Area are presented in Figs. C-11 through 
C-19 and the Husid plots for all the motions considered (except those at Emeryville) are 
shown in Fig. C-20. 


The time at which h(t) = 5 percent and the time at which h(t) = 95 percent are listed in 
Table C-3 for each component included in Figs. C-11 through C-19. Also listed in Table 
C-3 are the values of effective duration as defined in Appendix B. The range in effective 
duration for the horizontal motions recorded at these soft soil sites is shown in Fig. C-21. 
Also shown in Fig. C-21 are the values of effective duration obtained for the rock sites 
included in Appendix B (see Fig. B-40) and the values obtained using the equation 
relating effective duration at rock sites and distance (Eq. B-2). As can be noted, Eq. B-2 
appears to provide reasonable best estimates of duration for either site condition. 
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Table C - 1 
Peak Accelerations, Velocities and Displacements of Horizontal Earthquake Ground 
Motions Recorded at Soft Soil Sites during the Loma Prieta Earthquake 


Peak Peak Peak 


nae Acceleration | Velocity magi 

Station Component g cm/sec 
Be may eo NWOT 
saw [0.227 
ne oe 0.283 
0257 
San Francisco Airport | 633 0332 
0.235 


Alameda Naval Station 76.0 0.209 
0.268 
Outer Harbor Wharf 76.4 N55E 0.271 


Emeryville 51.0 
ais) 
Treasure Island 816 12.2 


Richmond City Hall- | 85.0 
Paiigiat | [Se as a 
Eeciaper ey Tenia | 99.0 BW [a7 oa 


Note: peak accelerations listed in this table are "Volume II" accelerations. 
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Table C - 2 
Values of v/ a and ad / v’ of Horizontal Earthquake Ground Motions Recorded at Soft 
Soil Sites during the Loma Prieta Earthquake 


Distance v/ a 
Station km Component | cm/sec/g ad/v’ 


APEEL No.2 1638 
Cho tage 
Foster City #83 
eS SC ee a 


OD SaaS me ee 
ese oTIIS dale 160. | 
cy | HL Sc EC ae ace 

Cae a 
Sts ee 
EW ise] 1d 


aS pans NT LCT NT TON 
Treasure Island 51.6 
peme—er eet 
Richmond City Hall- | 850 


Larkspur Ferry 
Terminal 
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Table C - 3 
Effective Duration Determined from Husid Plots of Horizontal Earthquake Ground 
Motions Recorded at Soft Soil Sites during the Loma Prieta Earthquake 


Effective 
Distance Time at 5% | Time at 95% Duration 
Station km Component seconds seconds seconds 


APEEL No. 2 46.8 
iy ek cle 
Foster City _ seeoe 
per | a eee ae 
San Francisco Airport LEW | 86s) | 1905 | 
Mee ea a 
Alameda Naval Station | 760 
ee Se 
Outer Harbor What | 764 
ia Sas 
Emeryville 51.0 
igo oe aa 
Treasure Island 51.6 
| 


Richmond City Hall- | 85. 
| a er 
Larkspur Feny 
Tew | oe 
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Fig. C-1 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at APEEL No. 2 during the Loma Prieta Earthquake 
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Fig. C-2 Accelerograms and Spectral Ordinates of the Horizontal Components of 
Earthquake Ground Motions Recorded at Foster City during the Loma Prieta Earthquake 
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Fig. C-3 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at San Francisco Airport during the Loma Prieta Earthquake 
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Fig. C-4 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Alameda Naval Station during the Loma Prieta Earthquake 
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Fig. C-5 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Outer Harbor Wharf during the Loma Prieta Earthquake 
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Fig. C-6 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Emeryville during the Loma Prieta Earthquake 
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Fig. C-7 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Treasure Island during the Loma Prieta Earthquake 
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Fig. C-8 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Richmond during the Loma Prieta Earthquake 
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Fig. C-9 Accelerograms and Spectral Ordinates of the Horizontal Components of Earthquake 
Ground Motions Recorded at Larkspur during the Loma Prieta Earthquake 
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Fig. C-10 Normalized Spectra of Horizontal Motions Recorded at Soft 
Soil Sites in San Francisco Bay Area During the Loma Prieta Earthquake 
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Fig. C-11 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at APEEL No. 2 during the Loma Prieta Earthquake 
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Fig. C-12 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Foster City during the Loma Prieta Earthquake 
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Fig. C-13 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at San Francisco Airport during the Loma Prieta Earthquake 


O 5 70 15 20 25 30 
Time - sec 
0.30 
0.15 
0.00 
-0.15 
NS Component 
-0.30 
O 70 20 30 
Time - sec 


Cumulative Normalized [a(t) *a(t}] 


Acceleration, a(t) - g 


O 5 70 15 20 25 30 
Time - sec 
0.30 
0.175 
0.00 
MEE? EW Component 
-0.30 
(6) 70 20 30 
Time - sec 


Fig. C-14 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Alameda Naval Station during the Loma Prieta Earthquake 
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Fig. C-15 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Outer Harbor Wharf during the Loma Prieta Earthquake 
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Fig. C-16 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Emeryville during the Loma Prieta Earthquake 
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Fig. C-17 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Treasure Island during the Loma Prieta Earthquake 
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Fig. C-18 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Richmond during the Loma Prieta Earthquake 
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Fig. C-19 Accelerograms and Husid Plots of the Horizontal Components of Earthquake 
Ground Motions Recorded at Larkspur during the Loma Prieta Earthquake 
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Fig. C-21 Effective Duration for Horizontal Motions Recorded at 
Soft Soil Sites and at Rock Sites During the Loma Prieta Earthquake 
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APPENDIX D 


INPUT DATA FOR EQUIVALENT LINEAR 
ANALYSES AT THE TREASURE ISLAND SITE AND 
AT THE SAN FRANCISCO AIRPORT SITE BASED ON 
BEST ESTIMATE SHEAR WAVE VELOCITIES 


To Report on: 


ASSESSMENT OF SITE RESPONSE 
ANALYSIS PROCEDURES 


APPENDIX D 


INPUT DATA FOR EQUIVALENT LINEAR ANALYSES AT THE 
TREASURE ISLAND SITE AND AT THE SAN FRANCISCO AIRPORT 
SITE BASED ON BEST ESTIMATE SHEAR WAVE VELOCITIES 


D.I INTRODUCTION 


The input data used for the equivalent linear analyses at the Treasure Island Site are listed 
in Section D.2 and those at the San Francisco Site are provided in Section D.3 of this 
appendix. These input data were used with computer program SHAKE91; the listings 
provided in this appendix follow the format specified in the user’s manual for SHAKE91 
(Idriss and Sun, 1992). Note that each sublayer thickness is given in feet, total unit 
weight is given in kip per cubic foot and shear wave velocity in feet per second. 
Thickness and shear wave velocity can be converted to meters and to meters per second 
by dividing the values given below by 3.28. Total unit weights can be converted to 
kN/m? by multiplying the unit weights given below by 157.2. 


D.2 INPUT DATA FOR THE TREASURE ISLAND SITE 


Option 1 -- dynamic soil properties - (max is thirteen) : 
1 
3 
abst #1 modulus for clay (seed & sun 1989) upper range 
0.0001 0.0003 0.001 0.003 0.01 0.03 0.1 O53 
abe Sie 10. 
1.000 1.000 1.000 0-981: 0.941 0.847 0.656 0.438 
0.238 0.144 0.110 
ital damping for clay (Idriss 1990) - (about Lower Range from SI, 1970 for sand) 
0.0001 0.0003 0.001 0.003 0.01 0.03 0.1 0.3 
ake 3.16 10. 
0.24 0.42 0.8 1.4 29 Sree 9.8 T5i55 
Pak Palsy ¢ 28. 
11 #2 modulus for sand (seed & idriss 1970) - upper Range 
0.0001 0.0003 0.001 0.003 0.01 0703 0.1 Os 
abe Sie Oe 
1.000 1.000 0.990 0.960 0.850 0.640 0.370 0.180 
0.080 0.050 0.035 
aka damping for clay (Idriss 1990) - (about Lower Range from SI, 1970 for sand) 
0.0001 0.0003 0.001 0.003 0.01 0.03 0.1 hee! 
ab. on al) 
0.24 0.42 0.8 1.4 2.8 Siu 9.8 1525 
lee 25 28. 
8 #3 modulus reduction for rock (Schnabel et al, 1972) 
-0001 0.0003 0.001 0.003 0.01 0.03 OF Ven 
1.000 1.000 0.9875 09525 0.900 0.810 On725 07550 
5 DAMPING IN ROCK AVERAGE 9/4 
-0001 0.001 0.01 0.1 ake 
0.4 0.8 ake 30 4.6 
3 1 2 3 
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Option 
2 


Ww OI nu RPh WD PP 


PRPPPR 
hWNHO 


15 
16 
17 
18 
19 
20 
7Aah 
22 
23 
24 
aS 
26 
749 
28 
29 
30 
31 
32 
33 
34 
35 
36 
Option 
= 


2 


3 


-- Soil Profile 


OWRPRPEPEPPPP PPP PNNNNNEPRPP PP RPP RP HPP NNNNNNNDND O 


1900 2048 


Option 
4 
36 
Option 
5 
0 
Option 


4 


TI - Best Estimate Shear Wave Velocities 


3.00 -050 -120 

3.00 -050 -120 

4.00 -050 -120 

5.00 -050 -125 

5.00 -050 ae 3) 

5.00 -050 ao 

5.00 -050 -125 

5.00 -050 -125 

5.00 -050 -125 

5.00 -050 -100 

5.00 -050 - 100 

5.00 -050 -100 

5.00 -050 - 100 

5.00 -050 - 100 

5.00 -050 - 100 

5.00 -050 - 100 

5.00 -050 - 100 

5.00 -050 -100 

5.00 -050 - 100 

5.00 -050 - 100 

10.00 -050 - 130 
10.00 -050 -130 
10.00 -050 - 130 
10.00 -050 - 130 
10.00 -050 - 130 
10.00 -050 - 130 
10.00 -050 - 130 
10.00 -050 -130 
10.00 -050 -130 
10.00 -050 - 130 
10.00 -050 - 130 
15.00 -050 -130 
20.00 -050 -130 
25.00 -050 - 130 
25.00 -050 -130 
-010 -140 


input motion: 


.02 yb_90.acc (8£10. 
-0671 20. 3 8 


800. 
800. 
800. 
500. 
500. 
525. 
525. 
550. 
550. 
600. 
600. 
700. 
700. 
700. 
700. 
600. 
600. 
600. 
600. 
600. 
1000. 
1000. 
1000. 
1000. 
900. 
900. 
900. 
900. 
900. 
900. 
900. 
900. 
900. 
1250. 
1250. 
4000. 


6) 


sublayer for input motion {within (1) or outcropping (0): 


number of iterations & ratio of avg. strain 


O75 
sublayers for which accn time histories are 


3 4 5 6 if 8 9 10 11 
1 1 1 1 1 1 1 1 1 
0 0 0 0 0 0 0 0 0 


sublayers for which accn time histories are 


18 19 20 21 22 23 24 Pats) 26 
al 1 z 1 1 at al 1 al 
0 0 0 0 0 0 0 0 0 


to max strain 


computed & saved: 


computed & saved: 


27 28 29 
1 al 1 
0 0 0 


Option 6 -- sublayers for which accn time histories are computed & saved: 
6 
30 ope 32 33 34 36 
Bi af 1 ik a af 


0 0 0 0 0 0 
execution will stop when program encounters 0: 
0 


D.3 INPUT DATA FOR THE SAN FRANCISCO AIRPORT SITE 


Option 1 -- dynamic soil properties - (max is thirteen): 
1 
3 
11 #1 modulus for clay (seed & sun 1989) upper range 
0.0001 0.0003 0.001 0.003 0.01 0.03 0.1 0.3 
alee Bie 10. 
1.000 1.000 1.000 0.981 0.941 0.847 0.656 0.438 
0.238 0.144 0.110 
ala damping for clay (Idriss 1990) - (about Lower Range from SI, 1970 for sand) 
0.0001 0.0003 0.001 0.003 0.01 0.03 OL 023 
ale 3.16 10. 
0.24 0.42 0.8 1.4 2.8 5.1 9.8 S35 
PAL PMV 28. 
akal #2 modulus for sand (seed & idriss 1970) - upper Range 
0.0001 0.0003 0.001 0.003 0.01 0.03 0.1 OR3 
ale She AOK 
1.000 1.000 0.990 0.960 0.850 0.640 0.370 0.180 
0.080 0.050 0.035 
allal damping for clay (Idriss 1990) - (about Lower Range from SI, 1970 for sand) 
0.0001 0.0003 0.001 0.003 0.01 0.03 OL 0.3 
ab Si. 10. 
0.24 0.42 0.8 1.4 2.8 SL 9.8 L5)-5 
Dave PHS. 6 28. 
8 #3 modulus reduction for rock (Schnable et al, 1972) 
-0001 0.0003 0.001 0.003 0.01 0.03 0.1 L.0 
1.000 1.000 0.9875 0.9525 0.900 0.810 0.725 0.550 
5 DAMPING IN ROCK AVERAGE 9/4 
-0001 0.001 0.01 0.1 Ae 
0.4 0.8 Li 3.0 4.6 
3 al 2 3 
Option 2 -- Soil Profile 
2 
1 36 SFO - Best Estimate Shear Wave Velocities 
1 2 5.00 -050 -130 1260. 
2 al 5.00 -050 -105 290. 
3 1 5.00 -050 .105 290. 
4 il 5.00 -050 -105 290. 
5 2 6.00 -050 eb20 500. 
6 2 6.00 -050 20 500. 
7, 2 8.00 -050 -120 1080. 
8 2 8.00 -050 -120 1080. 
9 2 8.00 -050 #120 1080. 
10 1 6.00 .050 -105 500. 
11 ak 8.00 .050 2125 1215. 
12 ab 10.00 -050 225 1215. 
13 ak 10.00 -050 225 12 Se 
14 1 10.00 .050 a25 125% 
15 il 10.00 .050 a2 ARPS Ushi 
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16 1 10.00 -050 125 sty. 
aly) al 10.00 -050 -125 ig aR 
18 1 15.00 .-050 -130 1680. 
19 a 15.00 -050 -130 1680. 
20 all 15.00 -050 ra .0 1450. 
+5 i. 15.00 -050 oho 1450. 
22 al 20.00 -050 -130 1950. 
23 al 20.00 -050 oi30 1950. 
24 x 20.00 -050 - 130 1950. 
25 1 20.00 -050 -130 2040. 
26 al 20.00 -050 -130 2040. 
27 1 20.00 705.0 -130 2040. 
28 1 20.00 -050 -130 2040. 
29 1 20.00 -050 -130 2040. 
30 at 20.00 -050 -130 2040. 
ay 1 20.00 -050 130 2040. 
32 1 20.00 -050 -130 2040. 
33 1 20.00 -050 -130 2040. 
34 ab 20.00 -050 230 2040. 
35 1 20.00 .050 130 2040. 
36 3 -010 -140 4000. 
Option 3 -- input motion: 
3 
1900 2048 .02 ssf.acc (8£10.6) 
- 1047 20. 3 8 
Option 4 -- sublayer for input motion {within (1) or outcropping (0): 
4 
36 0 
Option 5 -- number of iterations & ratio of avg. strain to max strain 
5 
0 8 0.5 
Option 6 -- sublayers for which accn time histories are computed & saved: 
6 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 
1 al ail 1 1 1 1 al 1 al al al af 1 
1 0 0 0 0 0 0 0 0 0 0 0 0 0 
Option 6 -- sublayers for which accn time histories are computed & saved: 
6 


16 17 18 19 20 aL 22 23 24 25 26 27 28 29 
1 1 1 al al 1 1 1 1 1 1 al a 1 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 
Option 6 -- sublayers for which accn time histories are computed & saved: 
6 


30 31 32 33 34 35 36 
1 ah 1 1 1 ak 1 


0 0 0 0 0 0 0 
execution will stop when program encounters 0: 
0 
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